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l/-} ; Disease outbreaks, such as those of Severe Acute Respiratory Syndrome in 2003 and the 2009 

pandemic A(H1N1) influenza, have highlighted the potential for airborne transmission in in- 
door environments. Respirable pathogen-carrying droplets provide a vector for the spatial 
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spread of infection with droplet transport determined by diffusive and convective processes. 
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An epidemiological model describing the spatial dynamics of disease transmission is presented. 
The effects of an ambient airflow, as an infection control, are incorporated leading to a delay 
equation, with droplet density dependent on the infectious density at a previous time. It is 
found that small droplets (~ 0.4 ^m) generate a negligible infectious force due to the small 
viral load and the associated duration they require to transmit infection. In contrast, larger 
droplets (~ 4 fim) can lead to an infectious wave propagating through a fully susceptible 
population or a secondary infection outbreak for a localised susceptible population. Droplet 
diffusion is found to be an inefficient mode of droplet transport leading to minimal spatial 
spread of infection. A threshold air velocity is derived, above which disease transmission is 
impaired even when the basic reproduction number Rq exceeds unity. 
^ ' Keywords: respiratory droplets, influenza, ventilation 
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1. Introduction 



l|, however, the spatial 



The temporal dynamics of disease transmission is well documented 
spread of infection and the mechanisms driving it are less well understood. Spatial dynamics 
of respiratory diseases has received much attention in recent years following the rapid global 
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spread of both Severe Acute Respiratory Syndrome (SARS) 2j and pandemic influenza (2009) 
A(H1N1) [3]. Large-scale geographic models are typically implemented by superimposing a 
transportation network on local infection dynamics 4j, [5[. Smaller scale models often imple- 
ment reaction-diffusion equations to describe random movements within populations [f], Q]. 
Such models have been shown to exhibit traveling-wave solutions whose existence depends on 
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the basic reproduction number R t 



However, few models address the mode of disease 



transmission that underlies the spatial dynamics. 

For some respiratory infections, such as influenza, three modes of transmission have been 
identified: airborne, droplet and contact transmission [10]. All three modes result from the gen- 
eration of respiratory droplets by an infected person during an expiratory event (e.g. coughing, 
sneezing). Droplet and contact transmission require relatively close contact between the in- 
fected and susceptible individuals for efficient disease transmission. Therefore, spatial spread 
via these routes of transmission must be driven by human movement. In contrast, the air- 
borne spread of a pathogen may be attributed to the movement of both people and fine aerosol 
droplets suspended in the air and their associated airborne residence time and pathogen load. 
Furthermore, the airborne route is the primary mode of transmission for other respiratory 
diseases, such as tuberculosis, and its contribution to the spatial spread of infection is of 
paramount importance. 

The relative importance of the three modes of transmission is difficult to quantify; however, 
recent experience with SARS and influenza outbreaks has highlighted the potential for airborne 
transmission in indoor environments. Evidence of the airborne transmission of respiratory 

n 

infections has been documented in hospital ward settings [ll(, housing complexes [12] and 
on board commercial airliners [IS]. In particular, the importance of indoor airflows has been 
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141 ] . with large air flow rates resulting in lower disease transmission [15f| . However, 



established 



the threshold air velocity above which disease transmission is impaired is as yet unknown. A 
comprehensive knowledge and understanding of the mechanisms driving such transmission is 
vital for the implementation of adequate infection controls in indoor public environments such 
as schools, hospitals and long-term care facilities. 

Airborne transmission is mediated by fine aerosol droplets small enough to remain sus- 
pended in air for prolonged periods, and large enough to contain non-negligible pathogen load. 



The standard epidemiological models operate on the assumption that contact between suscep- 
tible and infected persons is necessary for disease transmission [161 ]. However, the pathogen- 
carrying droplets emitted by an infected individual during an expiratory event are the disease 
vector and the standard models should be adjusted to reflect this. Such models have been 

n 

used to describe the spatial spread of fungal spores over a vineyard [17( and to model disease 
spread following the point release of an infectious agent [181 ] . 

A zero-dimensional model for the temporal development of an epidemic driven by expira- 



tory droplets (in particular, respirable droplets) was developed in 



19( | . The model is built on 



the concept of an infectious cloud surrounding each infected individual, an idea also considered 
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201 ] . In this work we investigate the airborne spread of an infection through a closed spatial 
environment within which the human population is confined for a prolonged period of time. Fo r 
example, such outbreaks have been documented in public settings including prisons 



boarding schools 



231 ] . long-term care facilities 



24 ] and on board large cruise ships 
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2jj|. In 



the absence of intervention measures, such outbreaks can persist for many weeks. We extend 
the zero-dimensional model to include the spatial dynamics of airborne-droplet transmission 
(also known as aerosol transmission) and investigate how an ambient airflow can influence 
disease spread. Accordingly, we neglect transmission by droplet spray, a close-contact trans- 
mission mode that occurs via direct deposition onto a susceptible's mucous membranes, and 
by (physical) contact transmission. Furthermore, we consider transmission only by respirable 
expiratory droplets, droplets whose post-evaporation diameter is less than 10 [im, neglecting 
transmission by inspirable droplets, droplets whose post-evaporation diameter is between 10 
and 100 //m. We follow [27|] and take the droplet post-evaporation diameter to be half the 



pre-evaporation diameter. Inspirable droplets contribute to disease transmission by inhalation 
almost immediately after generation (e.g., during the first breath) as they are considerably 
larger than respirable droplets and they gravitationally settle very fast. As in the case of 
droplet spray, transmission by inspirable droplets occurs only at close contact. 

2. A one-dimensional spatial model for airborne transmission 

Consider a population of susceptible, infected and recovered individuals. Let S(x, t) be the 
density of susceptibles, I(x, t) the density of infected and R(x, t) the density of recovered indi- 
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viduals, with N(x, t) = S + 1 + R being the total population density and n(to) = f N(x, to)dx 
representing the total number of people in the spatial domain at any time to ^ 0. Hence- 
forth, all densities refer to spatial densities, unless otherwise noted, and airborne droplets refer 
to respirable droplets. It is assumed that infected individuals continuously generate a cloud 
of pathogen carrying aerosol droplets and we let D(x,t) be the (number) density of active 
droplets, which we define as droplets that are both airborne and have a nonzero pathogen 
load. 

The zero-dimensional model, which constitutes the basis of the one-spatial dimension 
model, is derived in Stilianakis and Drossinos A general one-dimensional evolution equa- 
tion to model the spatial spread of disease due to the continuous motion of people and droplets 
takes the form 

dCi dQi _ , . 

IT— 5T + *'<*•'>• 

for each species i 6 [S, I, R, D], where Ci(x,t) represents the density of species i, the flux 
Qi(x,t) is the rate at which species i passes the point x at time t and ^i(x,t) is a source 
term representing the creation (or destruction) of species i. The density of susceptibles will 
decrease through contact with pathogen-carrying (respirable airborne) droplets and subsequent 
infection. The density of infected people will increase accordingly. Furthermore, it will decrease 
at the rate that individuals recover m, where is the disease infectivity period. Thus, 

q, s = -^DS, * I = ^DS-mI, Vr = ihI, 

per droplet of diameter d. A detailed derivation of the 



where /3d is the transmission rate 
transmission terms is provided in 191 ] ■ The dynamics of the airborne droplets is determined 
by generation and annihilation processes. The droplet density at any point is proportional to 
the density of infected individuals and increases at the rate Kd that pathogen-loaded droplets 
are shed. Moreover, the active droplet density will decrease as droplets are removed through 
gravitational settling and inhalation (by the person who generated it or another population 
member) and via pathogen inactivation. Thus, 

= Kdl - a d D, 
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28(] and the 



where ad is the droplet removal rate and l/ad the droplet infectivity period. 

We assume that human movement can be modeled as a diffusive process [g, 
associated flux is given by Fick's law as 

dC 

Qj = -D p -^, for j = S,I,R, 

where the minus sign is interpreted as the tendency of people to move from high density areas to 
low density areas and D p is the diffusivity of the human population. For real world situations, 
movement of the general population might also be modeled by including a convective term, 
whereby human motion would be faster and in a specified direction. However, movement in 
such facilities as prisons and long-term care facilities is more restricted and sporadic and a 
diffusive flux better models such a scenario. In addition, restricting movement to a diffusive 
process is more convenient to investigate if droplets can drive the transmission process. 

The droplet flux Qd depends on environmental characteristics. Droplets are generated 
through expiratory events (e.g. coughing, sneezing) with an initial velocity. Average expira- 
tion velocities are 11.7 ms" 1 and 3.9 ms _1 for coughing and speaking respectively [29]. For 
respirable droplets, the droplet relaxation time, the time required to adjust the droplet ve- 
locity to a new condition of forces, is of the order ~ 10~ 7 — s. Therefore, the droplet 
velocity rapidly tends to the carrier- gas (air ) velocity. The exhaled air flow has been modelled 



as a continuous turbulent round jet 
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311 ] . For a steady-state turbulent jet, small droplets, 
which follow the instantaneous fluid streamlines, may reach relatively large distances from the 
source (more than 8 m after 100 s). However, it was argued in 32)], where the air flow within a 



calm room was calculated via a Computational Fluid Dynamics software, that droplets smaller 
than 30/xm diameter would disperse within the room without significant influence of gravity or 
inertia. Accordingly, we consider that the initial velocity of the exhaled respirable droplets is 
the underlying fluid velocity, namely either the ambient air velocity, or zero in its absence. For 
convenience, we henceforth refer to the ambient airflow as the ventilation. This airflow could 
be naturally induced (e.g. a draft through an open window) or be the result of an artificial 
indoor ventilation system. 

It follows that, in an enclosed space with no ventilation, airborne droplets will be trans- 
ported by molecular diffusion alone and Qd = —Dd dD/dx, where Dd is the molecular droplet 



diffusivity. Molecular (Brownian) diffusivity may be safely ignored for droplets larger than 
0.5/um diameter. We retain it in the formal derivation of the flux equations to allow the sim- 
ulation of nanodroplets, to allow for turbulent droplet diffusion in the presence of an external 
turbulent flow, and to render the numerical solution of the parabolic differential equations 
stable. In a ventilated environment airborne droplets will also be convected by the ambient 
airflow. The appropriate flux is Qd = vD — D d dD/dx, where D d is the turbulent diffusivity, 
and v the average, constant air ventilation velocity. The one-dimensional droplet flux is an 
idealised approximation. In reality, the air flow is unlikely to be at a constant unidirectional 
velocity, and ventilation may simply mix droplet and air particles around within the domain. 
A higher dimensional spatial model is required to describe such complex flow characteristics. 
However, the one-dimensional model can be a good approximation for environments where air 
is predominantly in one direction, for example a draft of air blowing through a room. 
Assuming the diffusion coefficients of both people and droplets are constant the model is 

dS DS 8 2 S 



81 n DS d 2 I 



dN_ FN 

dt ~ Dp dx 2 ' [6) 
dD r n dD d 2 D 

— = Kd I-a d D-v— + D d —2, (4) 

and the density of recovered individuals is obtained from R = N — {S + I). Equations (Q3H|) 
reduce to the model equations proposed in [2] for infectious disease transmission by respirable 
droplets without contact transmission with zero ventilation velocity and diffusion coefficients. 

We consider the spread of a pathogen in the domain ^ x ^ I. We assume that people are 
confined to the interval [0, 1] and prescribe zero-flux conditions at the boundaries (the subscript 
x denotes partial differentiation with respect to x), 

S x (0,t)=S x (l,t)=0, I x (0,t)=I x (l,t)=0, N x (0,t) = N x (l,t) = 0. 

The boundary conditions for droplets would be expected to depend on the type of ventilation 
present. We assume that droplets are removed from the system when they reach the end of 
the domain (physically this could be, for example, through a wall vent or an open window). 
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One option is then to place an artificial boundary at the end of the domain x = I so that 
droplets effectively travel through the boundary unhindered. This is referred to as a transparent 
boundary condition. For convection-diffusion problems such a condition would take the form 

m, 

D t (l,t) + vD x (l,t) = 0. 

However, for a small diffusion coefficient, the boundaries are sufficiently far from the region of 
interest that, for the timescales of interest, droplets never reach them by diffusive processes 
alone. Therefore, it is sufficient to assume that droplets can be transported out of the domain 
by convection alone and we set the diffusive flux at the boundaries to zero, 

D x (Q,t) =D x (l,t) = 0. 

In our simulations we compared the use of the transparent and diffusive boundary conditions 
and found that, as expected, solutions were identical. Therefore, we present only those per- 
formed with the zero diffusive flux condition. 



The prescribed boundary conditions prevent the entrance or exit of people. Therefore, 
people are free to move about the domain but cannot leave it. Consequently, the total number 
of people in the domain will be constant for all time 

n (to) = / N(x,to) dx = no, V to > 0. 



Jo 

Initial conditions must also be prescribed to distribute people throughout the domain at t = 0. 
Since the total number of people in the domain is constant for all time, no is determined from 
the initial distribution 



n 



[ N(x,0)dx= [ \s(x,0) +I(x,0)l dx, (5) 
Jo Jo L J 



where we assume the initial number of recovered (immune) people is identically zero, R(x, 0) = 
0. Initially, we assume that no droplets are present and are generated for t > by the 
infected population. We consider two possibilities for the initial distribution of the human 

n 

population. The total population can be uniformly distributed throughout the domain 6], 
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which we henceforth refer to as a (spatially) homogeneous initial condition. This corresponds 
to a constant initial total population density 



In this scenario, N will be constant for all time, with N(x,t) = uq/1. Alternatively, people 
can be randomly distributed throughout the domain yielding a (spatially) heterogeneous initial 
condition. The latter case presents an opportunity to investigate whether droplets generated 
by infected individuals at one point in the domain can infect susceptible individuals at other 
locations, i.e., whether infection can occur without direct physical contact of a susceptible with 
an infected individual. 

2.1. Scaling and non-dimensionalisation 

We scale all population densities with the uniform density hq/1. The droplet scale is chosen 
by balancing droplet generation and removal processes such that, the number of droplets is 
approximately proportional to the number of infected people and the constant of proportion- 
ality is the ratio of the droplet generation and removal rates. We scale x with the length of 
the domain and choose the droplet removal time as the characteristic time scale, since this 
represents the time the droplet is airborne and capable of causing infection. Therefore, we 
scale 



AWioO^O) = — . 



S,I,N 




x ~ I 



1 

t r . 

ad 



(6) 



The dimensionless equations are 



dS_ 

dt 

m 

dt 
ON 
~dt 
8D_ 

~dt 




dx ^ d dx 2 ' 



(10) 



(7) 



(8) 



(9) 



where 



Ro — 



Pd^d 




V 



Vd, P = 



D d , P 



ad^i 



OLdV 



a d l 2 ' 
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The dimensionless parameters and the characteristic time scales of the model are summarized 
in Table [TJ The parameter Rq is the basic reproduction number, which describes the spread of 
disease through a completely susceptible population in the initial stages of an outbreak. A is 
the ratio of the droplet removal time scale r r to the disease infectivity time Tj and represents the 
fraction of the total disease infectivity period for which a droplet is capable of causing infection. 
The dimensionless coefficient of the convection term v = t t /t c represents the ratio of the 
droplet-removal time to the convection time. Similarly, the dimensionless number rja = T r /r,i 
represents the ratio of the removal time to the diffusion time. Boundary and initial conditions 
in dimensionless form are given by 

S x (0, t) = S x (l,t) = 0, S(x, 0) = S (x), 

4(0, t) = 4(1, t) = 0, I{x, 0) = J (a;), 

N x (0, t) = N x (l, t) = 0, N(x, 0) = S Q (x) + J (x), 

D x (0,t)=D x (l,t)=0, D(x,0)=0, 

where the functions So(x) and Jo(x) are prescribed initial population-density distributions. 
The dimensionless form of the initial population density, Eq. ([5]), is 



l f-i i 

dx = 1. 



/ N(x,0)dx=l S(x, 0) + I(x, 0) 
Jo Jo L 

For the homogeneous case, with constant initial population density Nh omo (x, 0) = 1, this 
implies 

S(x,0)+I(x,Q) = 1. 
3. Model parameters for an influenza outbreak 

We apply the one-dimensional model to numerically study the spatial and temporal dynam- 
ics of a model for an influenza epidemic. Most of the required parameters are taken from 
the parameters specific to the spatial dynamics and droplet size are discussed in the following, 
and parameter values used in the numerical simulations are summarized in Table [2j 

The effect of droplet size on disease spread is investigated by choosing two characteristic 
respirable-droplet sizes of post-evaporative diameters d\ = 4 /im and cfe = 0.4 fim. As men- 
tioned earlier, the post-evaporative diameter is taken to be half the emitted pre-evaporative 
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diameter. Droplet generation rates Kd are based on the number of pathogen loaded droplets 
emitted during a cough. Generations rates per cough are taken as = 160 day -1 



27J and 



Kd 2 = 240 day -1 [34]], the latter based on a cough volume of 400 cm 3 . The daily generation 
rates are obtained by considering a 200- fold increase for a sneeze [271 ] . and a total of 11 sneezes 
and 360 coughs per day [35]. 

In the model presented here, the infectious agent is not the droplet but the pathogens it 
carries. Therefore, the transmission rate per droplet /3d will depend the transmission rate per 
pathogen, /3d = [3 p qdNp, where qd is the probability of deposition in the human respiratory tract 
and Np(d) is the number of pathogens in a droplet of diameter d. In turn, the transmission 
rate per pathogen f3 p is determined from the contact rate Cd of a susceptible with a droplet and 
the probability pd that such a contact will result in successful transmission (3 p = CdPd- In order 
to derive the contact rate with a droplet we assume that each infected person is surrounded by 
a droplet cloud with volume V c i. It is further assumed that a susceptible individual comes in 
contact with a droplet through breathing during an encounter with this droplet cloud. If the 
average breathing rate is B and r ct is a characteristic time of breathing during the encounter 
then the contact rate Cd can be expressed as q = c^r rf , where c is the average number of 
total contacts a susceptible individual has per unit time. The transmission rate per droplet is 
thus 

Pd = c—T ct p d qdNp, 

and the number of pathogens per droplet can be determined by N® = VdP p , where Vd is the 
volume of the (spherical) pre-evaporative droplet and p p is the pathogen concentration of the 
lung fluid. Using the parameters of Table [2l the transmission rate per pathogen is calculated 
as 0.028 per day and the transmission rate per droplet thus evaluates to (3d 1 = 2.45 x 10~ 5 per 
day and j3d 2 = 5.57 x 10~ 9 per day. 

The droplet removal rate is determined by three distinct processes: gravitational settling 
of the droplet, inactivation of the pathogen load (which effectively removes the droplet) and 
the inhalation of the droplet by population members. Following Stilianakis and Drossinos 
the removal rate can be expressed as 

a d = d + Hp + (1 + CTa) y-Qd, (11) 
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where Od is the gravitational settling rate of the droplet and n p is the inactivation rate of 
airborne pathogens. Gravitational settling rates are size dependent and we take 6d 1 = 28.8 per 
day and 9d 2 = 0.39 per day [36|]. This implies that, in the absence of other removal processes, 
a droplet with diameter d\ = 4 ^m will settle under gravity in a tranquil environment in 
approximately 50 minutes. Decreasing the droplet diameter by an order of magnitude results 
in the droplet remaining airborne for approximately 2.56 days. The pathogen inactivation rate 
Hp, assumed to be independent of size, is taken to be 8.64 per day [33] and a droplet is effectively 
removed through pathogen inactivation in under 3 hours. Therefore, pathogen inactivation is a 
crucial process for the removal of smaller droplets that could, theoretically, take days to settle 
under gravity alone. Taking both gravitational settling and pathogen inactivation processes 
into account, droplets remain airborne and infectious for 38.5 minutes (d\) and 2.66 hours (d^)- 
The inclusion of removal through inhalation, the last term on the right-hand-side of Eq. (|lip . 
has negligible influence on these times (removal times of 35 mins or 2.46 hrs for d = d± or 
d = d2 respectively) and we approximate the droplet removal rate by = 9d + Hp- 

Droplet diffusivity will depend on the presence or absence of ventilation. In an unventilated 
environment, with v ~ m s~ 1 the molecular diffusivity of droplets is calculated as Dd 1 = 
6.2 x 10 -12 m 2 s -1 and Dd 2 = 8.34 x 10 -11 m 2 s _1 [36]. For an air-conditioned environment 



32] and 



with standard wall-mounted air-conditioners a typical airflow velocity is v = 0.2 m s _1 
we take this to be constant throughout the entire domain. Under such conditions the flow will 
invariably be turbulent and the droplet diffusivity will be several orders of magnitude larger. 
We estimate the turbulent diffusivity of both droplet classes to be D l d ur = 10~ 3 m 2 s _1 . 

Human diffusivity, motion, can be crudely estimated from D p ~ x 2 /t, where x is a charac- 



teristic distance traveled in a time t. For the spatial spread of a disease through a 



geographically 
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38, 



3&J. We 



open population typical distances traveled per hour are in the range 18 — 42 m 
expect that movement would be more restricted in a closed environment, for example a prison 
or long-term care facility, and we estimate that people diffuse approximately 10 meters per 
hour yielding D p ~ 10~ 5 m 2 s _1 . Intuitively, we estimate the spatial size of such environments 
to be of the order 10 3 meters and and fix our domain length at I = 2000 m. 
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4. Droplet dynamics 

4-1. Characteristic time scales 

The time scale of disease transmission n = 0C4/ ' (5d^d is determined from droplet properties 
and is thus dependent on droplet size. We estimate Tt(d\) ~ 3.72 days and r^cfe) ~ 2640 
days. Clearly, the time required by the smaller droplet to transmit disease will result in 
minimal disease transmission for the duration of the infectivity period (approximately 5 days 
for influenza). 

Droplet dynamics are greatly influenced by their (post-evaporative) diameter as this de- 
termines the residence time in air of an individual droplet, its infectivity properties and its 
pathogen load. The dynamics are described by three different timescales (Table Q]), two of 
which, T r (droplet removal time scale) and (droplet diffusion time scale), are dependent 
on droplet size. The convective timescale t c represents the time it takes for the ambient air- 
flow to carry the droplet the entire length of the domain. If I is of the order of 10 1 — 10 3 
m then typical convection times are r c = O(10 2 — 10 4 ) seconds and, clearly, convection is a 
fast process relative to the disease dynamics with r c <C 7*. Intuitively, this implies that lower 
ventilation velocities, which increase the convection time t c , result in greater disease transmis- 
sion. The droplet-removal (droplet infectivity) time scale is the inverse of the size-dependent 
droplet removal rate with T r {d\) ~ 38.46 minutes and T r {d,2) ~ 2.66 hours. Thus, since 
T r (di) < T r {d,2) we have that < z^ 2 and convection effects, for a constant air- ventilation 
velocity, are stronger for smaller droplets since they remain airborne longer and can be carried 
further from their point of origin. Furthermore, v is inversely proportional to the length I 
and convection will thus exert greater influence on droplet dynamics over shorter domains as 
droplets will be rapidly transported the entire length of the domain. 

The droplet diffusion time scales depend on the presence or absence of an ambient airflow. 
For an unventilated environment, typical Brownian diffusion time scales will be r&{d\) ~ 10 9 — 
10 13 days and r^cfe) ~ 10 s — 10 12 days for I ~ 10 1 — 10 3 m. In a ventilated environment, 
we estimate the time scale of turbulent diffusion to be T^ r ~ 1 — 10 4 days. Thus, for the 
short term spread of infection, since droplets are only airborne for 0(1 hour) and r r <C 7^, 
diffusion is an insufficient mechanism to transport droplets throughout the domain with the 
largest influence occurring at very short lengths and with non-zero airflow. However, over such 
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lengths contact and droplet transmission are more likely to drive disease transmission. 

From the preceding discussion on droplet timescales it is evident that, in a ventilated 
environment, droplets are rapidly transported out of the domain and disease transmission is 
comparatively too slow to result in substantial infection. For larger domains convection times 
can be significantly increased and the transmission of infection will become more efficient (e.g. 
r c w 3.5 days for I = 60,000 m). We estimate Rq = 1.34 and A = 0.005 for d\ = 4 /im and 
Rq = 0.0019 and A = 0.022 for d% = 0.4 /jm. The lower value of Rq for d = d<i is the result 
of the small viral load of the droplet which determines the magnitude of the transmission rate 
per droplet /3d, Tabled Thus, the larger generation rates associated with smaller droplets do 
not result in greater transmission. 

4-2. Droplet delay equation 

For a ventilated environment we can neglect droplet diffusion and Eq. fjXOj) becomes a first 
order equation. We consider the solution for the particular initial condition D(x,0) = 0. The 
characteristic curves satisfy 



which yields 



dt _ 1 
dx v ' 



t= ( X X *\ (12) 
v 



for any x = x* at t = 0. Along each characteristic curve D will satisfy 

dD 1 1 . 

_ + -J> = -1(3 

This first order linear equation can be easily solved for D{x) to obtain 

D(x) = — ( eh(Z,£=p) dC - / eh(£,i=p-) d£ 
v \Jo Jo 

The parameter be eliminated using Eq. (|12p yielding 

D(x,t) = - [ X e-^J^t-^de. 

v Jx-ut 

This solution is only valid far from the boundaries, for vt < x < 1, however, it provides an 
intuitive idea of the droplet dynamics. The role of droplets in disease transmission is thus to 
introduce a delay into the classic SIR system. Disease transmission at time t depends on the 
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density of infected individuals at a previous time t — ■ The integral also describes how 
the convection of droplets influences the spatial spread of the disease. It describes the force 
of infection that individuals at position £ exert on those at x and as \x — £| — > oo the droplet 
density vanishes and no infection will occur. Thus, the transmission risk decays with distance 
from the source of infection. A corresponding equation can be derived for the diffusion-only 
scenario showing similar behavior and its derivation is outlined in Appendix A. 

4-3. Stability of a homogeneous population density 

We derive stability criteria for the homogeneous case where the total population density 
is uniformly distributed throughout the domain and N = 1. For the basic disease- free state 
(S, I, D) = (1, 0, 0), the evolution of small perturbations (S, I, D) are governed by the linearised 
equations 

St = —XRoD + 7] p S xx , 
It = XR D - XI + r] p I xx , 
D t = I - D - ut) x + r] d D xx . 

We consider perturbations proportional to e tkx+ult , where u and k are the frequency and 
wavenumber, respectively. On substitution into the linearised system we obtain a quadratic 
equation for u 

u 2 + u{ivk + M + L) - XRq + ML + ivkM = 0, 
where, for convenience, we define 

M = X + T] p k 2 , L = 1 + r) d k 2 . 



This has solution 



vk (M + L) 

ijj(k) = —i 

v ' 2 2 



/ 4A.Ro - 4ML - v 2 k 2 vk(L - M) 
l±i H ttt tto 1- 2i- 



(13) 



(M + L) 2 (L + M) 2 

The disease-free state will be unstable when the real part ?R.{uj(k)} > and disturbances with 
wavenumber k will grow. Defining the square root in Eq. f)13j) as p + iq, with p > and both 
p and q real, we can write 
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Figure 1: 



where for instability we require p > 1. Now, equating real and imaginary parts yields 



21.2 



o n 4A.Ro - 4ML - v z k 

p — q = 1 + - 



vk(L - M) 



(M + L) 2 ' * p(L + M) 2 ' 

Solving for q in terms of p we obtain 

2 _ u 2 k 2 {L-M) 2 _ (M + L) 2 + 4A^p - 4ML - u 2 k 2 
{ - P '~ P p 2 {L + M) 4 (M + L) 2 ' 

Since F(p) is a monotonic increasing function of p the condition for instability can be expressed 
as F(jp) > F(l) which yields 



ML 



2 1.2 



1 + 



v k 



(14) 



(M + L) 2 _ 

We calculate r\ v = 5.77 x 10 -9 and % = 3.57 x 10~ 15 and, since % <C rj p <C 1, we can 
approximate Eq. (fl~4"|) by 



i?n > 1 + V 1 



k 2 



(A + l) 2 ' 

Clearly, Rq > 1 is a necessary condition for instability and it follows that the infection will 
always die out for d = cfe, regardless of whether ventilation is present or not, since i?o = 
0.0019 -C 1. In an unventilated environment, u = 0, the larger airborne droplets will be 
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transported by diffusion alone and the uniform state is unstable at all wavenumbers when 
Ro > 1 and the infection will spread. In a ventilated environment, v > 0, disturbances with a 
wavenumber satisfying 

k < ^^V«0-1 = fccrif (15) 

will be unstable. This result shows that the basic homogeneous disease-free state is stable 
to short wavelength perturbations (k — > oo) and only becomes unstable at long wavelengths 
(A; — > 0). The stability diagram is sketched in Figure [TJ Clearly, Rq must exceed unity for the 
onset of instability and, for fixed Rq > 1, the critical wavenumber decreases with increasing v, 
which physically corresponds to increasing the ventilation velocity and thereby reducing the 
time it takes for droplets to traverse the domain. Therefore, increasing v has a stabilizing 
effect, by reducing the number of unstable modes that will be amplified. 

For a velocity of v = 0.2 m s , the critical wavenumber is /c cr ;t = 2.54 and perturbations 
with wavelength satisfying A > A cr j t = ^~ ~ 2.47 will be amplified. Significantly, all per- 
missible wavelengths (A < 1) for our domain x £ (0, 1) are stable and therefore an arbitrary 
perturbation to the disease-free state will decay. Reducing the ventilation velocity decreases 
A cr it and unstable modes become permissible within the domain and any arbitrary perturba- 
tion will grow in time. The analysis indicates that, for a fixed size domain, increasing the 
ventilation velocity can mitigate the effects of a respiratory disease. 

5. Results 

We examined the spatial dynamics of the model for droplet size d\ = 4 //m considering both 
homogeneous and heterogeneous initial population-density distributions. The system ([7|)- (fTU]) 
was solved using a fully implicit finite difference scheme implemented in Matlab®. 

5.1. An unventilated environment 

In an unventilated environment, with u = 0, the spatial spread of disease will solely depend 
on the diffusion of people or droplets. Droplet diffusion is a slow process, for example a 4 ^m 
droplet will diffuse approximately 5.7 mm while airborne. In the same time a person could 
diffuse over 7 meters. Therefore, the spatial spread of disease will be driven by the diffusion of 
people as droplets are essentially stationary relative to the human population. Mathematically, 
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this observation follows immediately from the inequality rjd <C rj p . We first consider the disease 
spread when a localized infected density is introduced to a uniformly distributed susceptible 
population. We find that, for small times (t < 5 ~ 3 hours), droplet density rapidly increases 
until D(x,t) w I(x,t) and subsequently droplets and infected dynamics are indistinguishable. 
This occurs once droplet generation and removal processes are essentially balanced. This indi- 
cates that there is no separation between droplets and infected individuals and close contact 
is required for efficient transmission. At larger times the localized infectious peak is observed 
to grow in amplitude, Figure E|a). An outbreak in a closed population could be expected 
to persist for many weeks. However, times in excess of t = 2000 ~ 53 days are unrealistic 
for an influenza epidemic. At such times the model predicts the formation of two infectious 
traveling pulses, propagating in opposite directions. Accordingly, the density of susceptibles 
evolves into a wave front slowly infiltrating the completely susceptible population ahead of 
the front, Figure [2jb). There are several reasons why this behavior is observed. Firstly, the 
idealized case of a homogeneous population, where the infected population is surrounded by a 
constant source of susceptibles is uncharacteristic of a true human population distribution. In 
the absence of susceptible individuals the density of infected individuals would rapidly decay, 
Figure Secondly, once an outbreak is identified within the closed population, interventions 
such as isolation, quarantine or treatment would be imposed to limit the impact of the out- 
break. Diffusion models used to describe the larger scale geographic spread of disease display 
such traveling wave behavior, with transportation networks simulated using larger diffusion 
coefficients and timescales of the order of years being considered [{], Q] • 

Disease spread resulting from a heterogeneous initial distribution, where localized densities 
of infected and susceptible individuals are placed in different regions of the domain, is con- 
sidered in Figure [3l As before, the droplet density rapidly approaches the density of infected 
individuals and the two densities are thereafter indistinguishable from each other. We find 
that, even with very close contact between the two groups, diffusion is too slow a process to 
effectively transmit infection. In the absence of a susceptible population, and since r] p <C 1, the 
density of infectious individuals at any point x = xq will decay via I(t) = Io(xo)e _A< — > as 
t — > oo. The lifetime of the infected population is thus j which (dimensionally) corresponds to 
the disease infectivity time scale Tj of 5 days. The infected population recovers before diffusion 
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Figure 2: 



has time to effectively transmit the infection, r, <C r p . To conclude, in the absence of an 
ambient airflow, diffusion is an inefficient mechanism for the spatial spread of disease which 
will, presumably, be driven by contact and droplet transmission during close contact events. 

5.2. A ventilated environment 

In a ventilated environment, with v > 0, the infection dynamics are determined by the 
droplet size, the ambient air velocity and the length of the spatial domain. For a fixed droplet 
size and domain length, the wavenumbers that satisfy Eq. (115j) are determined solely by the 
ventilation velocity v. As discussed in Section 14.31 when v = 0.2 m s _1 , all permissible 
wavenumbers are stable. The evolution of an arbitrary perturbation under such conditions is 
shown in Figure [U Droplets generated by the infectious population are rapidly transported 
out of the domain causing minimal infection as they travel through the susceptible population, 
since the convective timescale is much less than that required for transmission r c <C Tt, and 
the infectious curve decays. 

If we rewrite (|15p . we find that the basic-state will be stable to an arbitrary perturbation 
provided the ventilation velocity satisfies 

a d l(\ + 1) r- - 

v > VRq - 1 = v ciit . 

ZTT 

For parameter values listed in Table [2] we find that f cr ;t = 0.08 m s -1 and ventilation velocities 
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of this magnitude and above are sufficient to prevent the spatial spread by the airborne route 
for a uniformly distributed population. However, we emphasize that this is an approximate 
value obtained for the relatively small diffusion of people and neglecting the possibility of con- 
tact and/or droplet transmission. In addition, the true air velocity will depend on the distance 
from the ventilation source, with larger velocities closer to the source. The unstable evolution 
of an arbitrary perturbation when v = 0.01 m s~ 1 is shown in Figure[5j An infectious pulse can 
be seen to propagate through the susceptible population in the direction of positive air flow. 
Droplet and infected densities are qualitatively similar and only slightly out of phase with each 
other, with droplets propagating ahead of the infectious pulse. Dynamically this implies that, 
for small air velocities, the dynamics of droplets and infectious individuals are closely coupled 
and only separate from each other at sufficiently high velocities. 

Finally, we consider the spatial transmission of infection resulting from a heterogeneous ini- 
tial population distribution, where infected and susceptible populations are placed at separate 
locations in the domain, Figure [6) For small times, t < 10 ~ 6 hours, the infected population 
generate droplets which are subsequently transported in the direction of positive air flow. A 
secondary curve for infectious individuals forms when droplets encounter the susceptible popu- 
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lation, t = 100 « 2.5 days. At later times these secondary cases produce droplets of their own 
and a secondary peak in the droplet curve can be observed, t = 500 ~ 13 days. The amplitude 
of the secondary curve for infectious individuals is obviously influenced by the initial distance 
between the infectious and susceptible populations. At large times, t = 1000 « 26 days, it 
is clear that the amplitudes of both curves for infectious individuals decay. This behaviour 
results from the decay of the initial infectious population through recovery and the subsequent 
reduction in droplets being convected towards the susceptible region. This demonstrates that 
it is possible for infection to occur in the absence of an infected population and purely as the 
result of the airborne non-diffusive transport of aerosol droplets. 

6. Discussion 

In this work we developed a model for the spatial spread of an airborne infection driven 
by pathogen-carrying droplets. Two droplet diameter values, d\ = 4 /im and c?2 = 0.4 /xm, 
were considered that are representative of experimentally determined droplet size distributions. 
Droplet dynamics are governed by generation and removal processes, the latter being domi- 
nated by gravitational settling and pathogen inactivation. Inactivation is a particularly crucial 
removal process for smaller droplets that could remain airborne for days. Smaller droplets 
are found to be a weak disease vector. Their relatively large generation rates do not result 
in greater transmission due to the small viral load and the associated duration required to 
transmit infection. Transmission is, thus, dominated by the larger droplets. 

A delay equation was derived for the droplet density as a function of the infected population 
density. The role of droplets in disease transmission is to introduce a delay into the system, with 
disease transmission at a given time dependent on the number of infected individuals present 
at a previous time. The equation also highlights how droplet-driven transmission decays with 
distance from the source of infection. 

The relative importance of diffusive and convective processes in the spatial spread of infec- 
tion was investigated. Two initial population distributions were considered: spatially homo- 
geneous or heterogeneous. In both cases droplet diffusion is shown to be a slow process with 
disease spread, in an unventilated environment, driven by human movement. This result follows 
from the observation that the time require for droplets to diffuse is significantly greater than 
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that for humans, t p <C r^. In the homogeneous scenario and for long time scales, the model 
displays the classic infectious wave propagating through a susceptible population, following 
an initial transient state until a balance is achieved between droplet generation and removal 
processes. In contrast, human diffusion is an insufficient mechanism to transmit disease in 
the heterogeneous scenario, with the infected population recovering before encountering sus- 
ceptible individuals, Tj <C t^. However, the inclusion of ventilation effects and the subsequent 
transport of droplets from the source of infection can result in a secondary outbreak if suscepti- 
ble individuals are encountered. This signifies that infection is possible without direct contact 
between susceptible and infected individuals. Furthermore, it was shown that increasing the 
velocity above a critical value can impair disease transmission in a homogeneously distributed 
population as droplets will be rapidly transported out of the domain causing minimal infection 
since the time required for transmission is large relative to the convective time scale r c <C rj. 

The use of ventilation to prevent disease transmission is well-accepted by both society 
and science. However, little work has been done on implementing control strategies based 
on this knowledge. Identifying optimum air velocities for indoor environments could mitigate 
transmission and reduce the disease burden in health-care facilities, schools and other densely 
populated locations. The model presented here indicates that, in a fixed size environment, the 
distribution of sick/healthy individuals and the ambient air velocity are the primary factors to 
consider when analysing such an intervention. 

Appendix A - Solution of droplet equation with diffusion 

In the absence of ventilation effects the droplet equation (jlOp takes the form 



and, for convenience, we consider the solution on an infinite domain as droplet diffusion is a 
sufficiently slow process that boundaries do not interfere with the solutions, Figures [2] and [3j 
We let D(x,0) = 4>(x) denote the initial droplet density and assume a solution of the form 



3D 
~dt 



D + I + Vd 



8 2 D 




dx 2 



D(x,t) = H(t)G{x,t) 
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which yields 



dG G(x,t)dH „. . I(x,t) d 2 G , rr/ . , 
+ ^7^— = -G(x,t) + 4-^ + % ^- T , with H(t)^0. 



dt H(t) dt 
We choose H(t) to satisfy 



~dt 



-H, 



H{t) Ia dx 2 



H{t) = e 



and the function G(x,t) then satisfies the inhomogeneous diffusion equation 

dG r , . a 2 G 

where f(x,t) = -fey an d with initial condition G(x,0) = 0(x). Following Kevorkian 
solution for G(x, i) will be of the form 



4QJ, the 



G(x, t) = u(x, t) + v(x, t), 



where u and v are the solutions of the simplified diffusion equations 



U t = VdUxx + f(x,t), 

u(x,0) = 0, 



Vt = rjd,v xx , 
v(x, 0) = (f>{x) 



These equations can easily be solved to obtain 



u{x, t) 
v(x, t) 



ex P ( - ^ 71 X 7a ) f( X 'i dz'cft', 



o 7-oo y / 4^(t - t') V 4 %(i - * 



i exp 



and the solution of (lA.lj) is then 



D(x,t) =e" 



t /"OO 



1 iap r-. (a :^^g^ 
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time scale (s) dimensionless parameter 

Tj = 7- disease infectivity 

7> = disease transmission Rn = = Ii 

1 Pd K d u a d m r t 

T r = — droplet removal (droplet infectivity) A = ^ = ^ 

t c = 1 - convective v = ^ = % 

Td, P = ifc diffusion (droplet, person) r) d , p = |g = ^- 



Table 1: Characteristic time scales and derived dimensionless parameters. 
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parameter 

Hi infection recovery rate 

c contact rate 

p p pathogen concentration in the lung fluid 

B breathing rate 

Vd personal-cloud volume of an infected person 

p d infection probability by an inhaled pathogen 

T ct characteristic breathing (contact) time 

P P transmission rate per inhaled pathogen 

v air velocity 

D d ur turbulent diffusivity of droplet 

Dp diffusivity of people 

Hp pathogen inactivation rate 

I domain length 

Parameters dependent on droplet size 

d droplet diameter (post-evaporation) 

Vd pre-evaporation (spherical) droplet volume 

Np initial number of pathogens per droplet 

q d inhaled droplet deposition probability 

Pd transmission rate per inhaled droplet 

Kd droplet production rate 

9 d gravitational settling rate 

a d droplet removal rate 

Dd molecular diffusivity of droplet 



value 

0.2 per day 
13 per day 

3.71 x 10 6 pathogens cm" 

24 m 3 d^ 1 

8 m 3 

0.052 

20 min 

0.028 per day 

0.2 m s" 1 

10~ 3 m 2 s" 1 

10~ 5 m 2 s" 1 

8.64 per day 

2000 m 



di=4/im 

d2=0.4/im 

V dl = 2.68 x lO" 10 

V d2 = 2.68 x 10~ 13 

JVp(di) = 9.95 x 10 

N p (d2) = 9.95 x 10 



cm 

-4 



5 per day 
9 per day 



0.88 
0.2 

= 2.45 x 10" 
= 5.57 x 10" 
= 4.1 x 10 5 per day 
= 6.14 x 10 5 per day 
28.80 per day 
0.39 per day 
= 37.44 per day 
= 9.03 per day 



Qd 2 

Pd 2 
Kd! 
Kd 2 
dd! 
0d 2 
Ctd! 

et d2 

D dl = 6.2 x 10~ J 
D d2 = 8.34 x 10' 



Table 2: Parameter values 
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Figure captions 

Figure 1: Stability diagram for a homogeneously distributed population. 

Figure 2: Model dynamics of airborne-influenza transmission driven by the diffusion of peo- 
ple, with droplet diameter d\ = 4//m. Initial conditions are Iq{x) = O.Ole^ 2 ^ -0 ' 5 ) 2 , 
Sq(x) = 1 — Iq(x) with wavenumber k = 50. (a): Density of infected individuals, (b): 
Density of susceptible individuals. All variables are dimensionless and scaled following 

Figure 3: Model dynamics of influenza transmission driven by the diffusion of people, with 
droplet diameter d\ = 4/im. The graph shows the densities of susceptible (solid line) and 
infected (dashed line) individuals with initial condition Iq(x) = ^L e - k2 ( x -° A ) 2 ) Sq(x) = 
25 e -k 2 (x-o.5) 2 anc | waV enumber k = 50. All variables are dimensionless and scaled 
following ©. 

Figure 4: Model dynamics of an influenza outbreak driven by the convection of droplets, with 
respirable droplet diameter d\ = 4/zm and ventilation velocity v = 0.2 m s . Initial 
conditions are I Q (x) = 0.01e _fe2 ^ _0 - 2 ^ 2 , S (x) = 1 - I (x) with k = 30. (a): Density of 
infected individuals (b): Density of droplets. All variables are dimensionless and scaled 
following dSJ). 

Figure 5: Model dynamics of an influenza outbreak driven by the convection of droplets, with 
respirable droplet diameter d\ = 4/xm and ventilation velocity v = 0.01 m s . Initial 
conditions are Iq(x) = O.Ole - ^^ -0 ' 2 ) 2 , Sq(x) = 1 — Iq{x) with k = 30. (a): Densities of 
infected individuals (solid line) and droplets (dashed line), (b): Density of susceptible 
individuals. All variables are dimensionless and scaled following ©. 

Figure 6: Model dynamics of an influenza outbreak driven by the convection of droplets, 
with respirable droplet diameter d\ = 4^m and ventilation velocity v = 0.2 m s . 
The graphs show the density of susceptibles (solid line) on the right vertical axis. The 
densities of infected individuals (dashed line) and droplets (dotted line) are shown on the 
left vertical axis. Initial conditions are Jo (a?) = ^ e - k2 ( x -° A ) 2 } Sq(x) = ^L e - k2 ( x -°s) 2 
and with k = 50. All variables are dimensionless and scaled following ([6]). 
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